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0.1 Introduction 

Entanglement is a pure quantum property that may appear when we have a composite 
object. It features the existence of a very special kind of correlations, which cannot 
occur for product states or mixtures thereof. Those correlations, sometimes referred 
to as quantum correlations, give rise to a great variety of phenomena and form the 
basis of many applications in quantum information science ([1} . This is why, during the 
last fifteen years, a large theoretical effort has been devoted to define and characterize 
this intriguing property. More recently, some of the ideas developed in this context 
have been used to get a novel perspective into the many-body quantum systems that 
appear in Nature. In particular, the role of entanglement has been analyzed in many 
papers, and its implications have been used to introduce novel ways of describing 
such systems. In this paper we review some of the basic concepts coming from the 
theory of entanglement, their applications to many-body systems, and some of the 
new theoretical methods that have come up in this context. 

Entanglement appears whenever we have two or more quantum objects. It high- 
lights the appearance of certain kind of correlations which cannot appear in classical 
theories. By object we may mean a particle, a bosonic or fermionic mode, etc. Strictly 
speaking, the Hilbert space corresponding to the whole system must decompose as a 
tensor product of several Hilbert space, one for each and every object. Entanglement 
will depend on how we understand this decomposition and thus we must specify with 
respect to what (particles, modes, etc) we refer when we discuss entanglement prop- 
erties of systems. Entanglement can also be defined for both pure or mixed states. In 
the first case the definition of entanglement possesses no difficulty and can be very 
naturally understood in terms of standard correlations. In the second, however, the 
definition is subtle as correlations do not directly imply entanglement (in fact, cor- 
relations are everywhere in the classical world!). For two objects, entanglement can 
be easily defined, qualified and quantified. For many objects, however, this is not the 
case. New possibilities (and correlations) may appear, which cannot be reduced to the 
case of two objects. In fact, as of today, the theory of entanglement has not been fully 
developed: there exist many open questions and some definitions need to be sharp- 
ened. In the first two sections of this chapter we attempt to give a rather superficial 
introduction to this topic, in order to prepare the reader for the following sections. We 
will first consider bipartite systems (ie, two objects) and later on multipartite ones. We 
will introduce separately entanglement for pure and mixed states. We will highlight 
two quantities: the entropy of entanglement and the quantum mutual information. 
The first one measures the bipartite entanglement for pure states, whereas the second 
one measures correlations for both pure and mixed states. A thorough review of the 
theory of entanglement can be found in 

In the next section we apply the ideas previously introduced to many-body quan- 
tum systems as they appear in typical physical scenarios. In particular, wc will consider 
spin lattices with short-range interactions and in thermal equilibrium. We will first 
show that certain entanglement measures display special features whenever a phase 
transition occurs, and explain that this is due to the simple fact that those measures 
basically display correlations, which are known to change abruptly under those circum- 
stances. A thorough review of the behavior of entanglement under phase transitions 
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can be found in ([3]). Then we will consider area laws in those systems. For that, wc 
will use the quantities previously defined, namely the entropy of entanglement and 
the quantum mutual information. When we select a region of the lattice and consider 
the existing entanglement and correlations between that region and its complement 
one finds that they scale with the number of spins lying at the boundary of such a re- 
gion. We will explain where this peculiar property comes from, and derive it for Gibbs 
states. There are many reviews of the area law and its application, like for example 
(HE]) 

In the last Section we will use the intuition developed through the area law in 
order to introduce an eflicient way of describing many-body quantum systems. It 
is based on a construction in which each spin in the lattice is replaced by several 
auxiliary spins, which are maximally entangled to their nearest neighbors. A map is 
then applied to each lattice site which transforms the auxiliary spins into the original 
ones, giving rise to the quantum many-body state. The resulting states are known 
as projected entangled-pair states (PEPS), and provide such an effective description. 
We will highlight some of their properties, and present some explicit examples in ID, 
where PEPS reduce to the so-called matrix product states (MPS), which play a central 
role in renormalization procedures and algorithms. Reviews on PEPS, MPS and other 
tensor network states can be found in ([HI [7]) 

In this chapter, we will not make any reference to the physical system we are 
dealing with. This is very standard in quantum information, whereby the theory is 
valid for all physical quantum systems with some reachable quantum levels. Thus, 
many of the ideas reviewed here may well apply to atoms in optical lattices, spins in 
magnetic materials, or electrons in solids. 

0.2 Entanglement in many-body systems: Pure States 

We start out considering the simplest case, namely the entanglement present in a 
many-body quantum state, 5', whenever we have a pure state. This is the case, for 
instance, at zero temperature (if there is no degeneracy), or in most applications in 
quantum information science. Unfortunately, it is an idealized case since zero temper- 
ature cannot be reached in practice and, in most of the systems, interaction with the 
environment leads to decoherence which reflects itself in the state becoming mixed. 
Nevertheless, pure states are much easier to deal with than mixed states, and already 
display many of the features entanglement is characterized for. In this section we will 
concentrate in this case and consider first the situation in which we only have two 
systems. Later on, we will analyze the more general case of multipartite systems. 

0.2.1 Bipartite systems 

We consider here two systems, A and B. We denote by Ha and Hb the corresponding 
Hilbert spaces, and by and orthonormal basis in Hx, where n = 1, 2, . . . , djc, 

with dx = dim(iJx)- Most of the time we will concentrate on the simplest systems, 
qubits, where cIa — ds ~ 2. In that case, in order to keep the standard notation, we 
will take as a basis {|0)x, |l)x}- Unless we state it differently, we will always work 
with qubits. 
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The Hilbcrt space corresponding to the whole system, H, is the tensor product of 
Ha and Hb, which we write H = Ha <E} Hb- An orthonormal basis in that space is 
{|n)A (8) \'m)B}- To simphfy the notation, we will typically omit the symbol "(8)", and 
the subindices A, B whenever there is no possible confusion. For instance, any state 
for two qubits can be written as 

1 1 

I*) = C„,™|n,m), \Cn,rrf = l. (0.1) 

n,m— 1 n.m—Q 

We will also omit the limits in the sum whenever it is obvious. 

We will consider observables for each of the systems, which will be represented 
as operators acting on the corresponding spaces. For instance, ® cr^ denotes an 
operator ai acting on A and a2 on B. As before, we will omit the symbol for the tensor 
product; additionally, will stand for af^l-^, where 1 is the identity operator. Pauli 
operators acting on qubits will often appear. They are defined as follows: 

a, = |0)(1| + |1)(0|, 

a,=*(|0)(l|-|l)(0|), 

a, = |l)(l|-|0)(0|=-ia,(Tj„ 

We say that 5* S Ha<E)Hb is a product state if there exist two vectors Lpi G Ha and 
(p2 G Hb such that |^) = \^i) a®\^2) b- Otherwise we say that ^ is an entangled state. 
Examples of product states are those forming the orhonrmal basis \n,m). Examples 
of entangled state are the so-called Bell states 

|$±) = i=(|0,0)±|l,l)), (0.2) 

|vl/±) = i=(|0,l)±|l,0». (0.3) 

Entanglement and correlations. The difference between entangled and product states 
is that the first ones give rise to correlations. If we have a product state l^*) = \ipi,ip2), 
then the expectation value (*|crf (8i cril"^) factorizes into {tpi\ai\ipi) {ip2\cri\ip2) , and 
thus the results of measurements in both systems will be uncorrclated. For any entan- 
gled state, on the contrary, there always exist observables in A and B for which the 
expectation value does not factorize and thus for which the results of measurements 
will be correlated. For example, considering the Pauli operator along the direction 
defined by a unit vector (in the xz-plane), n, as afi ~ n^cTx + n-zaz, we have 

{<i>+\a^<E><jg\<i>+)=n-m, (0.4) 

whereas ($+|cr4|$+) = 0. Thus, whenever we measure the same Pauli operators in 
A and B (ie, when n and rri are parallel), the results are random but completely 
correlated (ie the same outcome in A and B). 

In order to highlight the power of the correlations arising in entangled states, one 
can consider the following game ([S]). Two people, Alice and Bob, after meeting in 
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order to discuss their strategy, are isolated in two rooms and given each of them a 
different number, x,y ~ 0,1. Their goal is to output another number, a,b ~ 0,1, 
respectively, such that a (B b = xy (where © denotes addition modulo 2). That is, if 
(x, y) = (0, 0), (0, 1) or (1, 0), they must output a = b, whereas for (x, y) — (1, 1) they 
must give a different output. It is very easy to show that since Alice does not know the 
number y that Bob will receive (and viceversa), she does not know which a she has to 
output. For instance, if Alice is given x = I, Bob could have y = or ?/ = 1, in which 
case she would have to give a different value of a. But since she does not know y, there 
is no way they can always guess. The best strategy can be easily shown to give the 
right answer with probability P = 0.75. Now, if Alice and Bob during the discussion of 
the strategy share an entangled state $+, then they can guess with a probability P = 
0.853... > 0.75. The strategy is: (i) for Alice to measure on her qubit the observable (Tfi^ , 
whenever she is given x, where uq = (0,1), ni = (1,0); (ii) for Bob, the observable 
cTjjij^ whenever he is given y, where rrio ~ ( 1 / \/2, 1 /a/2) and mi = (1/a/2, — l/v^). 
It is very simple using ()0.4|) to determine that the probability of guessing is indeed 
(1 + V2) / {2^/2) . Thus, in a world where entangled states exist, one can do things that 
are otherwise impossible, like to play this game with a higher probability. Actually, this 
is the idea behind quantum information science, where by using quantum states and 
correlations one can perform tasks (in the context of cryptography or computation) 
that are classically impossible. Furthermore, if one performs this game experimentally 
(for instance, using photons) one obtains a result which is incompatible with local 
realistic theories (ie, theories where the outcome of measurements do not depend on 
what is being measured somewhere else, and where the properties we measure are 
already well defined previous to the act of measuring). This is the essence of Bell's 
Theorem ([9]), which states that local realistic theories are incompatible with quantum 
mechanics. Note that the Bell states are the ones that gives a highest probability of 
wining the game, and thus they give rise to the maximal quantum correlations. They 
are usually referred to as maximally entangled states. 

Schmidt decomposition. In order to analyze the entanglement in bipartite systems 
it is useful to introduce the Schmidt decomposition (SD). Given the state ^E* (|0.ip . it 
is always possible to find an orthonormal basis {!««)} in Ha and {|vrn)} in Hb, such 
that we can write 



The reason is that the matrix C can always be written as C = UDV, where U and 
V are isometries and D is diagonal with positive elements. This way of writing C is 
called singular decomposition, and it is valid for any matrix (|10p. When (Ia = ds the 
square matrices U,V are unitary {UU'^ = U'^U = 1). Otherwise they are rectangular 
but still fulfill the required properties such that we can always find the Schmieedt 
decomposition. Note that U and V can be found by diagonalizing CC^ and C^C, 
respectively. The matrix D can be found by taking the square root of the resulting 
diagonal matrix, which coincides in both cases. The diagonal elements of D are the dk 
that appear in (j0.5p and are called Schmidt coefficients. 

Once equipped with the SD, we can easily figure out wheater a state is entangled 
or not. In case one of the Schmidt coefficients is one and the rest are zero, we have 




(0.5) 



k<.dA -ds 
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a product state. Otherwise, our state is entangled. The first statement is obvious, 
whereas the second follows immediately when one tries to write the state as a product 
state. 

The SD is also very useful to determine the reduced density operators for subsys- 
tems A and B alone. Using (|0.5|) wc find 

PA^tTB{mm = Y,di\uk){uk\ 

k 

PB = ivAmm = Y.<^l\vk){vk\ 

k 

where by trx we mean the trace with respect to system X. The reduced density opera- 
tors so obtained are automatically diagonalized, and we observe that their eigenvalues 
arc nothing but the square of the Schmidt coefhcients. 
Let us consider the following example: 

I*) = cos(6i)|0, 0) + sin(6')|l, 1), (0.6) 

for 6* G [0, 7r/4]. This is already written in the SD form. The eigenvalues of the reduced 
density operator are cos^lO) and s\t?{0). For = we have a product state, whereas 
for = 7r/4 we have a Bell state, which is the one that gives the largest quantum 
correlations, as explained in the context of the game. In parallel, the reduced density 
operators get more and more mixed as one increases 9 from to 7r/4 (note that 
the purity of a mixed state, p, is related to the distribution of its eigenvalues when 
considered as probabilities). 

Entropy of entanglement. The previous example indicates that entanglement is re- 
lated to the mixedness of the reduced density operators. In fact, this relation sug- 
gests that we can introduce a measure of entanglement by using any measure of 
mixedness of a state. A very natural way of the latter is the von Neumann entropy 
S{p) — — tr[plog2(jo)] (in the context of quantum information, one defines the loga- 
rithm in base 2). Thus, we define the entropy of entanglement of a state (jlip . VP, 

E{-^) = -tT[pA\og^{pA)] = -tr[pB\og^{pB)] = -Y,dl\og^{dl). (0.7) 

k 

As we see, this quantity can be easily determined through the SD. Note that in this 
expression log(O) = by definition. For a product state, = 0, whereas the maxi- 
mum entanglement is E ^ log[min(dA) (^b)]; which is reached for the state for which 
all the dk are equal. Those are thus called maximally entangled states (even for d> 2). 

In order to give a physical interpretation to the entropy of entanglement, we will 
now define two key concepts in the context of quantum information: entanglement 
concentration and distillation ([T^ . But, before doing that, we make a parenthesis 
to introduce yet another basic feature in quantum information and which is needed 
in order to explain the afore mentioned concepts, namely generalized measurements. 
The postulates of Quantum Mechanics specify what occurs when we perform a so- 
called filtering measurement of an observable. One just has to consider the associate 
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Fig. 0.1 (a) Entanglement distillation. Out of many copies of a weakly entangled states, by 
LOCC we obtain less copies of a maximally entangled state (b) Entanglement dilution. The 
reverse process. 

operator, say O, and its spectral decomposition. Let us denote by Oi its eigenvalues, and 
by Hi the projector onto the corresponding eigenspaces. In case Oi is not degenerate, 
Ilj = \4>i){(f)i\, where 0\(j)i) = Oi\(f)i). Note that X^^n = 1. Given a, in general, mixed 
state represented by a density operator p, the probability of obtaining the outcome Oi 
is Pi = tr(ni/9), and the state after the measurement is pi = HipWi/ Pi. Thus, a filtering 
measurement is characterized by a set of operators Hi which are positive (since they 
are projectors, their eigenvalues are zero or one) and add up to the identity operator. In 
fact, there is a more general scenario representing a measurement. We may considering 
bringing a measuring apparatus, letting interact with our system, and then reading 
off the apparatus. According to the laws of Quantum Mechanics, we will describe the 
measurement as follows. First, the initial state of the system and the apparatus will 
he p® |^)(^|, where \ A) denotes the (pure) sate of the latter. The interaction can be 
described in terms of a unitary operator, [/, so that the state becomes U {p®\A) {A\)U'^ . 
Denoting by = \A^j){Af^\ the projector operators defining the measurement on the 
apparatus, we have that the probability, P^, of obtaining the outcome labeled by /i, 
and the state of our system after the measurement, p^, can be written as 

AnpA], 

p^=tr(^^p4), = (0.8) 
where A^ = {A^\U\0) is an operator acting on our system. Note that 



Y,AlA,^{A\U 



C/t|A)=l, (0.9) 



where we have used the fact that \A^) is an orthonormal basis. Thus, a generalized 
measurement is described in terms of a set of operators, A^, fulfilling (|0.9|) . In fact, it 
can be easily shown that for any set of operators fulfilling that condition, there exists a 
measurement scheme (ie states \A) and | A^)) which renders ()0.8p ((Ij. Here, for the sake 
of simplicity we have restricted ourselves to pure states of the measurement apparatus 
and non-degenerate measurements. The extension to other cases is straightforward. 
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Once we have defined generalized measurements, we can introduce the concept of 
entanglement distillation [Fig. lO.lf a')]. Let us first consider the state (|0.6p . and let 
us assume that our goal is to create a maximally entangled state (ie, a state in the 
same form but with 6 ~ 7r/4) by acting locally on each of the particles. In order 
to do that, we can try to apply a generalized measurement to the first particle. We 
choose Aq ^ tan(6')|0)(0| + |1)(1| and Ai = (1 - AJAq)^/^. Those two operators 
fulfill (|0.9|) and this ensures that there is a physical measurement scheme associated 
to them which can be easily determined. In case we measure and obtain the outcome 
associated to Aq we will achieve our goal. Otherwise, we will produce a product state 
instead. The probability of succeeding is po = 2sin^(0). Note that, if A and B are 
spatially separated and are hold by Alice and Bob, respectively, in order to know if 
the measurement has been successful, the outcome (ie a classical bit of information) has 
to be transmitted from Alice to Bob. One says that by local operations and classical 
communication (LOCC) one can distill a maximally entangled state out of the state 
^' with probability pq. One may wonder if there is another generalized measurement 
applied to A and B (individually) giving a higher probability of success. In fact, this 
is not the case since the generalized measurement we just chose is the optimal one. 
Now one can consider the case in which Alice and Bob possess two identical copies of 
the state Vl', and they try to obtain maximally entangled states by LOCC (in which 
joint measurements on both qubits of Alice, or both qubits of Bob, are authorized). 
In general, they may get as outcome a maximally entangled state in a Hilbert space 
of dimension d = 2,3,4. For instance, if they are completely successful, they will get 
two copies of a maximally entangled state which is equivalent to a single copy of such 
a state in a space of dimension d = 4. Or if they get a single copy, they will have 
d = 2. One can show that, indeed in the case of two copies the average entanglement 
is strictly larger than twice that for a single copy if one chooses the optimal strategy. 
Now we can consider what happens when we take n copies and allow for the optimal 
LOCC in order to optimize 



where pd is the probability they end up with an entangled state in a space of dimension 
d (for d = 1 they end up in a product state) . The logarithm is the right quantity such 
that n copies of a maximally entangled state (which corresponds to a dimension 2") 
exactly gives a factor n. It turns out that in the limit n — oo the result precisely 
coincides with the entropy of entanglement E{'^) (|0.7p . This occurs not only for qubits, 
but for any d-level systems. Thus, the entanglement entropy is nothing but the optimal 
averaged entanglement that we can distill out of by LOCC in the asymptotic limit 
where we have a large number of copies. 

One may consider the opposite process, called entanglement dilution [Fig. lO.lf b)]. 
Given n maximally entangled states (of qubits), and by applying the optimal LOCC, 
how many copies, m, of the state ^' we can obtain. The ratio D = n/m (in average) 
in the limit n — )■ oo turns out to coincide again with E{'^), giving again a physical 
meaning to the latter. In fact, in this limit, states \E'i,2 can be converted into each other 
with a yield E {"^f 2) / E {"i/ 1) . This implies that, at least when we consider an scenario 




(0.10) 
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where we dispose of many copies of a state and are allowed to perform LOCC, there 
is just a single measure of entanglement, namely the entropy of entanglement. 

We finish this section by mentioning other quantities that are usually employed 
to quantify entanglement. One is the concurrence ()14p . which is the square of the 
determinant of the matrix c in ()0.1|) and is nothing but the product of the Schmidt 
coefficients. Another one is the fidelity with a maximally entangled state, ie 



where the maximization is with respect to the unitary operators U and V. It measures 
in a sense how close we are to a maximally entangled state; U and V just correspond 
to a basis change. 

0.2.2 Multipartite systems 

Entanglement in multipartite systems becomes more complicated than in bipartite 
ones. First of all, one can have that certain objects are entangled to others, but not 
to all of them. Second, the quantification becomes more subtle since it is not known 
if a property like the inter-convertibility of states by distillation and dilution exists. 

We say that a state 5' of systems A, B , . . . , Z , is a product state if there exists 
\ipx) € Hx such that j^*) = ^x\'Px)x- Otherwise we say that we have an entangled 
state. Still, it may happen that some of the systems are disentangled. In order to 
characterize the entanglement we consider all possible partitions of the systems, and 
for each of them we apply the above definition. Thus, entanglement is characterized in 
terms of a partition which indicates which systems are entangled among themselves. 
For instance, for three parties we can have: (i) they are in a product state; (ii) only 
A and B are entangled; (iii) only A and C; (iv) only B and C; (v) all are entangled. 
These cases are mutually disjoint. An example of case (i) is the state |0, 0,0), of case 
(u) \^+)ab ® |0)c, and of case (iu) the states (US] [HI) 



(we have omitted the normalization). 

These two last examples illustrate the difficulty of quantifying entanglement in 
many-body quantum systems. It is not clear which of those states is "more" entan- 
gled. In the context of quantum information, it depends on the application we have in 
mind. For some, one of them is more useful, whereas for other, it is the other one. Fur- 
thermore, in this case it is not possible to convert the state \ W) into the state \GHZ) 
by LOCC in the asymptotic limit without loosing copies (equivalently, entanglement), 
and thus we cannot assign a quantity like the entanglement entropy to them ()17p . One 
could hope to be able to convert n copies of any state ^ into mi copies of \W) and 
7712 of \GHZ), and then back in the limit n — >■ cxj, in which case one could define two 
measures of entanglement, given by the ratios mi. 2/^1. However, this is also impossible. 
One may then try to include other representative states beyond those two and thus 
define more measures of entanglement. However, it is not even clear that a finite set 
of representative states exists. Thus, one cannot follow the procedure we reviewed in 



= max|($+|([/^®V^B)|*)|2, 



(0.11) 



\W) = |0,0,1) + |0,1,0) + |1,0,0) 
\GHZ) = |0,0,0) + |1,1,1) 



(0.12) 
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the previous subsection to assign a meaningful entanglement measure to multipartite 
states. 

What we can still do is to consider bipartite partitions, in which we consider the 
entropy of entanglement of two disjoint sets of subsystems. This, in fact, has interesting 
applications in the context of many-body physics at zero temperature and will be 
analyzed in more detail in section [0.41 

Another approach is to look at fidelities with respect to certain particular states. 
For instance, one can define similarly to (|0.1ip the fidelity with respect to a GHZ, a 
W state; or to products of Bell states. Yet another possibility is to define measures 
which give figures of merit in specific applications. For instance, in the context of 
quantum repeaters, one usually defines the localizahle entanglement (|18p as follows. 
We consider we perform measurements in all particles except for two and the goal 
is to obtain, in average, the maximum entanglement for those to particles. Let us 
consider four particles, A, B,C, D, and define the localizable entanglement of \I/ with 
respect to particles A and D. Assume we measure certain observables Ob,c in B and 
C, respectively, and denote by \oi) an eigenbasis of those operators. Then, the average 
entanglement we obtain in A and D if we perform the measurement will be 



where Pij is the probability we obtain the outcomes i and j in B and C, respectively, 
and the state of AD in such a case. This entanglement will, in general, depend on 
the observables we decide to measure. We thus define the localizable entanglement as 
the maximum with respect to all possible observables. As mentioned before, this has 
a specific meaning in the context of quantum repeaters whereby the goal is to obtain 
as much entanglement as possible between the first and the last node by measuring in 
the intermediate ones. 

0.3 Entanglement in many-body systems: Mixed States 

In this section we consider mixed states. Those are described by a density operator, 
p, fulfilling (0 = > (meaning that all eigenvalues are non- negative) , and tr(p) ~ 1 
(normalization condition). It can always be written as 



where Pi > fulfill pi = 1 and the normalized states 4'^ need not be orthogonal. 
The interpretation of (|0.14p is the following: the state p can be obtained by preparing 
the state with probability pi (and then "forgetting" which state has been prepared). 
If one of the pi = 1 we have a pure state, which obviously fulfills tr(p^) = 1; otherwise, 
we say that our state is mixed. Note that for mixed states there exist many ways of 
writing (j0.14p . That is, there exist other qi > and |$i) such that p can be decomposed 
in their terms. For instance. 




(0.13) 




(0.14) 



p-^m{o\ + \i){i\) = -{\+){+\ + \-){-\), 



(0.15) 
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where |±) = (|0) ± |l))/\/2. This simply means that the same state can be prepared 
in different ways, either by mixing |0) and |1), or by mixing |±). 

The probabihty of obtaining an outcome fi when performing a (generahzed) mea- 
surement is given by ()0.8p . Thus, the density operator contains all information we have 
about our system. This entails that we will not be able to distinguish by any means 
how we prepared a state (which decomposition we used) given that the statistics of 
any measurement do not depend on how we prepared the state. 

In this section we will consider the entanglement of a mixed state p. We will 
start out with the simplest case, that of two subsystems, and later on analyze the 
multipartite case. 

0.3.1 Bipartite systems 

We consider again two subsystems, A and B. We say that p represents a product state 
whenever we can find pa,b, operators acting on Ha,b such that p = pA® Pb- ^ state 
is separable ()19p if p can be written as a mixture of product state; that is, if 

p ^^Pi\ai,bi){a^,hi\ (0-16) 

i 

where pi > 0. Otherwise, we say that p represents (or is) an entangled state. 

The definition of entangled states is self-explanatory. A product state is obviously 
not entangled, since it does not yield correlations. A separable state cannot be en- 
tangled either, since it is a mixture of product state, which themselves are not. But 
separable states may contain correlations; and this is what differentiates the pure and 
mixed states. For example, the state 

p=i(|0,0)(0,0| + |l,l)(l,l|) (0.17) 

fulfills {erf af) = 1 whereas {crf'^) = 0. These correlations are, however, very 
trivial. If we had a classical system we could also have them. Only entangled states 
may display non-classical correlations. 

A subtle point is that a state may look entangled even though it is separable. Let 
us take, for instance 

a = i(|<D+)(cf>+| + |a>-)(<f>-|) (0.18) 

According to this formula, we can prepare a by mixing two maximally entangled states. 
However, there is another way of preparing the same state a which does not require 
using entangled states at all. This immediately follows from the fact that a = p (one 
just has to replace the definition of the Bell states in this formula), and thus according 
to ()0.17p one can prepare it by mixing two product states. The state tr is thus separable. 
This simple example illustrates the difficulty of finding out whether a state is entangled 
or not. We have to check all possible decompositions: only if none of them involves 
product states we will have an entangled state. Unfortunately, there exist infinitely 
many decompositions, so that this task is hopeless. Fortunately, in some special cases 
there are shortcuts which can give us the right answer with much less effort. 
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Entanglement witnesses. An entanglement witness ^IU\\ is an observable which de- 
tects (witnesses) the presence of entanglement. Given an operator W = W"^ we say 
that it is a witness if for all product states |a, &), {a,b\W\a,b) > 0, but W possesses 
negative eigenvalues. From the definition of separable state (|0.16p . it is clear that if 
ti{pW) < 0, then p must be entangled. Thus, a negative expectation value of a wit- 
ness indicates the presence of entanglement. Note, however, that the converse is not 
necessarily true: if a the expectation value of a witness is positive, this does not imply 
that the corresponding state is entangled. However, one can show that for any entan- 
gled state there always exists a witness that detects it. Sadly, there is no simple way 
of finding out such a witness, which makes the problem of detecting entanglement in 
mixed state rather non-trivial. 

Let us consider some examples. First, for two qubits, W = 2 — 5*, where S = 
cr^^® (erf -l-cr^)-l-cr^<8)(o-f — cr^) and the sigma's are Pauli operators, is an entanglement 
witness. This can be shown by noting that for a product state, | (5*) | < | (erf) -I- (cr^) | -|- 
|(crf ) — (cr^)l < 2, since (a) < 1. By choosing the sigma's as in the game in Subsection 
[0T2T1 we see that for p{p) = (1 - p)l/4 -t- p|<I>+)($+|, ti[p{p)W] = 2^/2p, and thus 
entanglement is detected for p > 1/ V2. 

For our second example, we choose continuous variable systems. In this case, Ha b 
are infinite dimensional [and isomorphic to L^(i?)]. We consider two canonical opera- 
tors, X and P, for A and another two, Y and Q, for B. They fulfill canonical commuta- 
tion relations [X,P] = [Y,Q] = i. We define W = 5-2 with S = {X -Y)^ + {P + Q)^. 
Foraproduct state, (S) = AX^ + AY^ + AP^ + AQ^ + {{X) - (Y))'^ + {{P) + (Q))^ > 2, 
as a consequence of Heisenberg uncertainty relation AX'^ + AP^ > 1 (and similarly 
for Y and Q). On the other hand, since [X — Y, P + Q] =0, there always exist states 
for which both quantities are as small as we want, so that W becomes negative. For 
instance, two mode squeezed states are detected by this witness. 

Partial transposition. Another way of detecting entanglement is via the partial trans- 
position (j2ip . Given p, we can always write it in terms of an orthonormal basis {|7i, m)}, 
as 

P=Y. P^,r,kAhk){],l\. (0.19) 

We define its partial transpose with respect to A in the basis {|n)}, 

= E Pw\j^k){i,l\. (0.20) 

i.j,k,l 

For separable states, is a valid density operator (hermitian and positive) as it 
can be directly shown from ()0.16|) . However, for entangled states this need not be the 
case. Thus, if p^"^ has any negative eigenvalue the p must necessarily be entangled. In 
particular, for two quibts it can be shown (|20p that for all entangled states, the partial 
transposed is no longer positive semidefinite, which provides us with a very powerful 
tool to detect entanglement. For higher dimensional systems, however, there typically 
exist entangled states where p"^^ is positive semidefinite. In those cases, the detection 
of entanglement can be very complex as one has to find an appropriate witness. 
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As an example, let us consider the state p{p) defined in Subsection (|0.3.ip . A neg- 
ative eigenvalue of p{p)'^^ appears whenever p > 1/3, and thus this state is entangled 
in such a ease and separable otherwise. 

Entanglement measures and mutual information. For mixed state, entanglement mea- 
sures similar to the entanglement entropy can be defined. However, they are much 
harder to evaluate. For instance, we can consider the distillation procedure as before 
but now with mixed states (j^^ . That is, we may try to distill out of n copies of a state 
p the maximal number m of states using LOCC. The ratio m/n in the limit n — >■ oo 
is called distillable entanglement, D{p). Analogously, we may define the process of en- 
tanglement dilution and define the entanglement cost, Edp). In general, D{p) < Ec{p) 
so that we cannot first distill and then get back the same state as before. Further- 
more, there exist very few examples where those quantities can be evaluated. Another 
measure that can be determined in practice (at least for qubits) is the entanglement 
of formation 

EF{p)^mmJ2p^Ei■9), (0.21) 

i 

where the minimization is done with respect to all decompositions of p [cf (|0.14l) ]. This 
quantity is related to the entanglement cost through Ep{p'^"')/n Ec{p) in the limit 
n — > oo. 

Another way of measuring the entanglement is through the fidelity with a max- 
imally entangled state, as in (|0.1ip but now with F{p) — niax($+|(J7 (8) V)p{U'^ (E) 

Finally, another way of measuring entanglement is using the definition of partial 
transposition. One defines the negativity ([23|) as N{p) = max (| jp^'* | |i — 1,0) where 
the 1-norm is given by the sum of the absolute values of the eigenvalues (one can also 
define in terms of the logarithm of such an expression). The negativity can be positive 
only if we have an entangled state; however there exist entangled states for which 
it is zero. Nevertheless, it possesses certain properties which makes it very useful to 
quantify entanglement in a simple way. 

Another quantity of interest in the context of quantum information is the quantum 
mutual information, I (A : i?) (jlj. This does not measure entanglement, but rather 
correlations. In fact, it is the finest measure of correlations in the sense that it detects 
them even when correlation functions do not. It is defined through 

I{A ■.B) = Sa + Sb- Sab (0.22) 

Here, Sx is the von Neuman entropy of p restricted to system X. Thus, Sab = 
— tr(plog2 /o), and Sa = — ti'(pA log2 p^), where pA = trB(p) is the reduced density 
operator of subsystem A. The mutual information has the following properties: 

I{A:B)>0; I{A : B) ^ Q ^ p = pA ® Pb, (0.23) 

i||p -pa®Pb\\1< I{A : B) < log2(d)||p -pa®Pb\\u (0.24) 

I{A : B) < I{{aA) : B) < I{A : B) + 2Sa. (0.25) 

The first indicates that it is only zero for product states, ie when there are no cor- 
relations. The second relates it to standard correlations: whenever there are strong 
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correlations, the mutual information is large; the converse is not entirely true, since 
the dimension of the Hilbert space, d, appears in the expression, which may be very 
large. Finally, the last line indicates that it decreases whenever we discard a subsys- 
tem (in this case a), but it cannot decrease by more than twice the entropy of such a 
system. The mutual information has the interpretation that is the amount of qubits 
we must erase in order to obtain a product state ([24| . The proof of all those properties 
(except for the second one ([T])) easily follows from the strong subadditivity of the von 
Neumann entropy: given three subsystems, X, Y, and Z, Sxy + Sxz > Sxyz + Sx- 
If we choose Y = A, Z = B and X a system which is disentangled, we immediately 
obtain the first property. By choosing a, X = A, Y = B and Z such that the whole 
state 4" is pure and the reduced state in AB is our state p, we obtain Sab > Sa — Sb, 
which can be used to prove the third property: taking Z = a instead we obtain the first 
inequality, whereas using SauB > Sab — Sa and Sa + Sa > SAa one readily obtains 
the second one. 

0.3.2 Multipartite systems 

The description of multipartite entanglement must still confront several challenges for 
mixed states. Whereas the definitions of product and separable states are straightfor- 
ward, now we have to consider again different partitions in order to characterize this 
eluding property. So, given a partition of all subsystems in disjoint sets, we say that 
this partition is entangled if we cannot write the state as a mixture of product states 
along each set (but still entangled within each set; compare with multipartite pure 
states, and bipartite mixed states). We end up with a table in which for each partition 
we state whether the state is entangled or not. In order to check the entanglement for 
each partition, we have to find the appropriate witness, whose definition follows very 
naturally that for bipartite systems. Such a table contains some redundancies since, 
for instance, if a tripartite state is separable with respect to the partition {A){B){C) it 
automatically is for any other partition. However, not all the implications which occur 
for pure stats concur for mixed ones. For instance, there can be tripartite states that 
are entangled with respect to the partitions {AB){C), {AC){B), but not for {BC)A 
([25l) . This simply means that the state may be prepared by acting together (ie letting 
them interact) on A and B, and independently on C (plus classical communication), 
or, alternatively, by acting together on A and C, and independently on B. But never 
if A is not allowed to interact with B or C. 

As measures of entanglement one can use the ones defined for mixed bipartite 
states properly extended to many subsystems, or those for multipartite pure states, 
like the localizable entanglement. However, typically one finds other measures which 
are more appropriate to describe specific experimental situations. In the following 
we briefly review the spin squeezing p6p . which plays an important role in certain 
precision measurements. 

Spin squeezing. We consider a set of N qubits, although it is simple to extend it to 
d- level systems. We denote by 5 = {Sx^Sy, Sz) the collective spin operators 

1 ^ 

n=l 
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These operators fulfill angular momentum commutation relations [Sa^ Sp] = 2i€a,p,-ySj 
(with e the antisymmetric tensor and a sum over 7 is understood). We define the spin 
squeezing parameter 

.-^1ML_ (0 27) 

It can be shown that whenever C < Ij the state must be entangled (P7j) . This quantity 
cannot be written in terms of an entanglement witness since it is not liner in the state. 
However, one can easily define in terms of a witness for two copies of the state. 

For separable states, ^ > 1. This inequality is saturated when all the qubits are 
in the same pure state along the XY plane in the Bloch sphere, for instance in |+). 
Actually, the Bloch sphere can help us to get an intuitive picture of the meaning of 
this quantity. There, given a state we represent the possible values of the measurement 
of the Sa by an ellipsoid centered at (S'), and with axis given by the corresponding 
variances I^Sx, ^Sy and A5'z. For the state |+)®^ we will have a circle in the YZ plane 
centered at the point {N/2, 0, 0) and with radius y/N/A, yielding ^ = 1. This quantity 
may decrease if we deform the circle into an ellipse stretched along the Y direction 
and squeezed along Z. In that case, we can decrease the value of ASz but keeping (Sx) 
practically constant (in reality it will decrease, but since it is close to N/2 3> 1 this 
will not affect much ^ as long as the deformation is small). Thus, ^ < 1 corresponds 
to states which are close to pure product states (ie they are close to the surface of 
the Bloch sphere in the XY plane) but for which the variance of is reduced with 
respect to product states (at the expense of increasing some other variance). In fact, if 
N is large and we are close to the Bloch sphere, we can use the formalism of Holstein- 
Primakoff (I28|) in order to reexpress the situation in terms of two canonical variables. 
Let us assume that we work with states for which (Sx) — N/2; then we can treat it as 
a c-number (ie neglect its quantum fluctuations which will be very small as compared 
to its expectation value), and define two operators X = 2Sy/N, and P = 2Sz/N. The 
commutation relation [X, P] , when acting on the states we consider can be replaced 
by [X, P] ~ i, and thus we are left with two canonical operators, which fulfill the 
uncertainty relation AX'^ + AP^ < 1. States fulfilling ^ < 1 are now states fulfilling 
AP^ < 1/2, ie the so-called squeezed states. 

Particle vs. mode entanglement. We close this section by illustrating something we 
have already commented in the introduction; namely, that the property of entan- 
glement depends on how we define our subsystems. Consider bosonic particles, each 
of them possessing two levels, |0) and We can first consider the Hilbert space 
H — Ha <E) Hb (E) . . ., where Hx defines the space for each particle. Obviously, since 
we have bosons, not all states in H are relevant, but just those which are symmetric 
under exchange of any pair of particles. On the other hand, we can consider Fork 
space Hp = ® Hi, where Hq contains states |n)o with n particles in state |0) (and 
similar for Hi). Again, if we have a fixed number of particles, not all the states will 
be relevant. In the first (second) case we say that we write the states in first (second) 
quantization. Let us consider two state written in both languages: 

IV'i) = ^(|0)a ®\1)b + \1)a ® \0)b) = |1)o ® |l)i, 
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IV2) = ^(|0)a + = ^(|l)o ® |0)i + |0)o ® |l)i), 

(0.28) 

The first state is entangled in first quantization but a product in second, whereas the 
second one is the other way round. In fact, in the second one there is only one particle, 
but still we can have entanglement if we express the states in terms of modes and thus 
divide our space into that of the first and that of the second mode. This illustrates 
that, in order to properly speak about entanglement, we have to specify what are the 
subsystems we are considering. 

0.4 Entanglement and area laws 

So far we have been dealing with many-body states and their entanglement properties 
regardless of whether they are in equilibrium or not. We did not specify any interac- 
tions, nor temperature, but just defined some of the tools that are required to study 
entanglement for any state, in equilibrium or not. In this Section we concentrate on 
many-body states in thermal equilibrium which interact with short range interactions 
in lattices. This corresponds to some of the most interesting situations in experiments 
with cold atoms and other condensed matter systems. The requirement that we con- 
centrate on lattices is in order to avoid some of the mathematical inconveniences of 
working in the continuum, although some of the ideas reviewed here may be extended 
to the continuum just by taking the lattice constant a to zero. 

Entanglement in many-body systems has been and is an active area of theoretical 
research in the last years. Of particular interest has been the behavior of such a 
property along phase transitions. We will first very briefly review some of those results, 
and then concentrate on a property, the area law, that seems to be fulfilled by all 
systems in equilibrium, as we will explain, and thus strongly characterizes many-body 
quantum states that appear in Nature. We will consider systems both at zero and finite 
temperature. In the first case, we will use the entanglement entropy and a bipartition 
of our system. In the latter, stronger statements can be made, where we will use the 
quantum mutual information. 

We consider a spin system on a lattice in d spatial dimensions. Some of the ideas can 
be extended to bosonic and fermionic systems; nevertheless, for the sake of simplicity 
we will use the spin language and comment on the possible extensions to fermions. 
Spins interact according to a Hamiltonian, H, which we will assume that only involves 
few spins (typically only two, ie we will have two-body interactions) close to each 
other. In other words, we will consider finite-range interactions only. We will also 
be interested in the thermodynamic limit where the number of spins, N ^ 00. We 
consider thermal equilibrium states at temperature T, ie, described by the density 
operator 

PT = (0.29) 

where Z = tr(e~^/-^) is the partition function (and we have used units with fc^ = 1). 
At zero temperature, the state p reduces to a projector onto the ground subspace, ie, 
that fulfiUing — E^y^i), where is the ground state energy. If the ground 
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state is not degenerate (which we will assume in many occasions), we will just use a 
pure state \^o) to denote the ground state. 

0.4.1 Entanglement and phase transitions 

In thermal equilibrium, the state (j0.29p will depend on external parameters driving 
the Hamiltonian, which we will denote with a single letter B, as well as on the temper- 
ature. B could be an external magnetic field, interaction strengths, etc. Thus, we can 
write p{B,T). As we change one of those external parameters the state may change 
abruptly its properties in the thermodynamic limit. Strictly speaking, some of the 
observablcs may depend non analytically on some of those parameters. At the point 
of non-analyticity we say that the systems undergoes a phase transition. Particularly 
interesting transitions may occur at T = when we change the external parameter B. 
In that case, we speak about a quantum phase transition, since quantum fluctuations 
are now responsible for the transition. 

If we study the reduced state of two particles, tr, as we undergo a (quantum) phase 
transition, we may experience an abrupt behavior in its entanglement (j29l I30p . The 
reason can be easily understood as follows. Imagine we use a measure of entanglement 
(like the ones introduced for bipartite mixed states, Section lD.3.ip . say E[a{B)\. Then, 
if a correlation function presents a non-analyticity with respect to J5, this will also 
show up as a non-analyticity in the entanglement. The reason is that a is an analytic 
function of all correlation functions, and thus whenever the latter is not analytic, then 
neither is a and thus E. For example (schematically), 

dE[a{B)] _ dE Sa 

dB ~SadB' ^"-^"^ 

can become discontinuous whenever da/dB is, or even when the measure of entan- 
glement is not analytic. In other words, if E is not analytic at some point then it is 
because of the presence of a quantum phase transition, or because the definition of the 
entanglement measure is not analytic at that point. Thus, in some sense, looking at 
the entanglement of two subsystems does not really give more information about the 
phase transition than looking at two-body correlations. One has to look at the whole 
many-body system. 

There exist two methods that look at the whole many-body state to study phase 
transitions. The first is the based on looking at quantities like F{Bo) = lim^^Q \ (\['(i3o + 
e)|^'(i?o))| properly re-scaled, since they display special features whenever there is a 
quantum phase transition (|31l |32)) . The other method consists of using the localizable 
entanglement introduced in Section 10.2.21 There, for instance, hidden orders can be 
identified (03]). 

0.4.2 Area laws at zero temperature 

In this Section we will consider T = and a pure ground state of H. We consider 
a connected region A of the lattice with a smooth boundary and the complementary 
region B and concentrate on the entropy of entanglement, Ea, between those two 
regions and its growth as we make region A larger and larger. According to ()0.7p . this 
is given by the von Neumann entropy of the reduced density operator corresponding to 
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Fig. 0.2 Only particles that are at a distance smaller than the correlation length, ^, are 
correlated (or entangled) to the particles in B. Thus, only they contribute to the entanglement 
entropy and therefore the area law, Ea < (.Nqa- 

region A, pA- In general, since the entropy is an extensive quantity, one would expect 
that it scales with the number of spins in region A. In fact, this is the case if we take 
a random state on the lattice (according to the appropriate measure). However, for 
ground states of Hamiltonians as we are considering here, this seems not to be the 
case. Instead, the entanglement scales with the number of particles at the border of 
region A, NgA- In a sense, the entanglement scales not with the volume of region A 
but with the area, and thus the name area law (|34p . 

Let us consider a simple example, the Majumdar-Ghosh (j35]) Hamiltonian 



where the spins are qubits and N is even. Each term in the sum can be written as a 
projector onto the subspace of two qubis with total spin 1 plus an irrelevant constant. 
Thus, the ground state is a dimerized state, where singlet sates l^*") are formed 
between nearest neighbors. In fact it is doubly degenerate: one can have dimers formed 
between spins 2i — 1 and 2i or between 2i and 2i + 1. If we chose a set of neighboring 
spins as region A, the entropy of entanglement is bounded by 2, since at most two 
singlets will contribute to it. As we see, this bound does not grow with the size of A, 
and the border of A has a fixed number of spins (two) which does not grow either. 

The intuition behind the area law is displayed in Fig. 10.21 Note that for continuous 
systems the entropy generally is unbounded since Nsa — > oo when we take the lattice 
constant a to zero. Nevertheless, by scaling it properly an area law also arises. 

In the last few years, the area law has been studied in a variety of systems. Here is a 
summary of some of the results (|4j[5]): In 1-dimensional chains, it has been proven that 
all gapped systems fulfill the area law (p6| . For gapless systems which can be described 
in terms of a conformal field theory, the area law is slightly violated: Ea < | log2(^), 
where £ ~ L/a is the number of spins in region A, and c the central charge of the con- 
formal theory (|571 155)) . Interestingly enough, this is a universal law since the coefficient 
in front of the logarithmic violation is a constant independent of the lattice constant. 
Many other models have been analyzed both analytically and numerically, and no vio- 
lation (beyond the logarithmic one) has been found, at least for all reasonable models 
(ie homogeneous, finite range interactions, etc). 




(0.31) 
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In higher dimensions, the situation is not so clear. It has been conjectured, and it 
seems to be corroborated by all examples studied so far, that for gapped systems the 
area law is fulfilled. For critical systems there is at most a logarithmic correction (p9l 
40|) , even though for some of them one can show that an area law is not violated (j4T|) . 
In 2-dimensions we can thus conjecture Ea — ci£ log £ + + C3 log™ £ + j+ where 
£ = L/a and L the length of the boundary, ci 7^ for certain critical systems (like free 
fermions), and C3 ^ for the rest. All coefficients Ci are now not universal, since they 
change if we change the lattice constant. However, 7 can indeed be universal and, in 
fact, it can be connected to the existence of topological properties (|42l |43)) . 

As an illustration, we calculate the entanglement entropy for a set of free bosons 
in 1-dimension. We consider R lattice sites, and a Hamiltonian for free Bosons with 
periodic boundary conditions 

R 

H ^ ~t^{alar+i + al^j^ar) (0.32) 



where are annihilation operators with standard commutation relations [a^, aj] = 6, 
and with afl+i = ai. The ground state is given by 



r,s 



N 

|*o) = ^|^> a, I I vac), (0.33) 



with |vac) the state with bosons. We separate now the R sites into the first L <^ R 
and the rest, and are interested in the mode entanglement fsee 10.3.2)) between these 
two subsystems. We thus separate the whole Hilbert space into the tensor product of 
the Fock space corresponding to the first L modes, and the rest. In this representation, 
we can write 



N 



1*0) = E 
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n/2 


R-L 




R 




R 



{N~n)/2 

\n,N~n). (0.34) 



From this expression it is simple to calculate the reduced density operator for the first 

n)(7i|. For A^ ^ 1 we have 

1*00 

El^- p{n)\ogMn)]dn, (0.35) 



subsystem, p = X]^!Lo?'('^)l"')("l- A^ » 1 we have 



with 

p(n) ~ _^^^(n-no f /2c^ (O.36) 

V 27rcr^ 

no = NL/R and = NL{R — L)/R?. Performing the integral, we obtain Ea = 
log2 (T+const, which in the limit N,R oo with N/R constant it scales like 1/2 logj L. 

0.4.3 Area laws at finite temperature 

For finite temperature, we can not longer use the entropy to derive an area law since: (i) 
it does not measure correlations anymore; (ii) it is an extensive quantity at any finite 
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T. However, by using the quantum mutual information introduced in Section [0.3. II we 
will be able to rigorously derive (|33]) an area law valid for arbitrary finite temperature 
T 7^ as long as we have finite-range interactions. Note that at zero temperature the 
mutual information reduces (up to a factor of two) to the entanglement entropy, and 
thus it provides us with a suitable generalization of such quantity. Nevertheless, our 
bound will diverge at zero temperature as it should be expected from the fact that 
for critical systems one obtains logarithmic corrections to the area law, as explained 
in the previous section. 

We will consider nearest neighbor interactions, which can always be achieved under 
finite range interactions by redefining a spin which accumulates those of region of size 
equal to the interaction length. We write the Hamiltonian as 

H = Y,h,., (0.37) 

<i.j 

where hij is a Hamiltonian acting on nearest neighbor spins, and which we assume 
to be bounded above < ft- by a constant. Given a region A as before and its 

complement we want to show that 

IiA:B)<NsA!^log2{e). (0.38) 

In order to show that, we will use the fact that the free energy F{p) = {H)p — 
T S (p) / log2{e) , where p is any valid density operator and S the von Neumann entropy, 
is minimized by the Gibbs state (j0.29p . In particular, F{pt) < F{pa(^Pb), where pa,b 
are the reduced density operators of regions A and B, respectively. Using the definition 
of the free energy we obtain 

liA ■.B)^Sa + Sb- Sab < ^^^^'^ -^^)p-®pb i^g^(^)_ 39) 

By noting that the expectation value of each hi_j with pT and pa ® Pb coincides 
whenever i and j are both in region A or _B, we arrive to the desired expression (j0.38p . 

We remark that this area law is generally valid for any dimension and short-range 
Hamiltonian. In fact, one can also admit sufficiently fast decaying terms in the proof. 
Note also that for T the bound diverges, at it should be according to our discussion 
of critical systems in the previous section. 

0.4.4 Area law and correlation length 

In Section (0321 we gave an intuitive picture of the area law in terms of the correlation 
length (see Fig. 10. 2p . In this section we will make this connection more precise using 
again the quantum mutual information ([44| • This result will be independent of the 
temperature. 

As we mentioned above, I{A : B) measures correlations between regions A and B. 
Thus, as those regions get further away, the mutual information must decay. Indeed, 
one of the properties of the mutual information ()0.25|) expressed exactly this fact. If 
we fix region A and separate more and more the region B (see Fig. 10. 3[) , the mutual 
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Fig. 0.3 (a) We choose two regions A and _B separated by a distance i. (b) The mutual 
information decay as a function of L for each choice of A. 

information will be reduced by a factor of two at some point, € [1, oo] . We define the 
correlation length ^/ as the maximum of all with respect to all regions A fulfilling 
the requirements we imposed above. Note that this quantity may be infinite, like for 
instance, in critical systems at zero temperature. Using (j0.25p and the fact that in 
between regions A and B there arc of the order of NqaL sites, each of them with a 
maximal entropy equal to log2(2s + 1) (s is the total spin), we have 

Il=o{A) < Il=o (A) + 2NgA^i log2(2s + 1) (0.40) 

Using now that Il=^j{A) < Il=o{A)/2 we finally arrive to the area law 

lL=oiA) < msA^i log2(2s + 1) (0.41) 

which gives a finite bond for finite correlation length. 

We remark that the upper bound diverges for critical systems, as it should be 
given the discussion in Section 10.4.21 Note that the converse may not be true: if the 
correlation length is infinite, we may still have an area law. 

0.4.5 Detecting the area law 

How can we measure the area law? This seems to be difficult since measuring entropies 
is not an easy task, given the fact that it does not correspond to any expectation value 
of a physical observable. However, we may aim at measuring Renyi entropies instead, 
which arc defined as Sa (p) = log2 (trp" )/(l — a) . In fact, for a integer this corresponds 
to measuring quantities like tr(p"), which can be carried out if several copies of the 
system are at our disposal. The von Neumann entropy can be then determined by 
analytic continuation. In any case, one can use the Renyi's entropies instead of the 
von Neumann entropy in most of our previous definitions, and expect similar behaviors. 

The quantity tr(p^) = tr[(p p)r®^], where T is the swap operator between each 
site and its copy partner, ie 

T = ^ |n,TO)(m,n|. (0.42) 

Thus, we just have to measure the observable T on each particle and its copy, multiply 
the results, and then average over many measurements. For qubits, T — Si ■ S2 up 
to a constant, and thus could be measured by letting the qubits interact for a while 
according to the Heisenberg Hamiltonian and then measure both of them. In a similar 
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way, one can measure the n-th Renyi entropy by using n copies. Note that one would 
expect that the outcome of the measurement would get exponentially small with A'', 
something which should happen for random states. However, in the presence of an 
area law, this decay with TV will be softened and, in 1 spatial dimension, may allow to 
carry out experiments with subsystems containing large number of sites. 

0.5 Tensor network states 

Many-body quantum systems are very hard to describe since the number of parameters 
required to specify a state scales exponentially with the number of lattice sites. The 
reason is that the Hilbert space of the whole system is the tensor product of those 
corresponding to each lattice site, and thus the dimension of that space displays the 
exponential scaling. In the previous section we have seen that there is a common 
property of thermal states in lattices systems with short range interactions, namely 
the area law. This fact can guide us to find out an efficient language to describe many- 
body quantum systems in which the number of parameters only scales polynomially 
with the number of sites. Under certain conditions, this is indeed possible in terms 
of so-called tensor network states. In this section we argue how the area law leads to 
such descriptions and briefly review some of them. 

0.5.1 Area law and Projected entangled-pair states 

Let us consider the Majumdar-Ghosh model (|0.31[) which clearly illustrates how the 
area law arises. Let us consider a dimerized ground states, as displayed in Fig. I0.4r a). 
There are two spins in each site, each of them maximally entangled with the neighbor- 
ing spins. The entanglement entropy of any region A with the rest, B, equals 2, since 
only the entangled states at the border of the region contribute to the entanglement, 
each of them with one unit. If the spin Hilbert space would have dimension D instead 
of 2, then we would have 2 log2(I?). Now, imagine we map each pairs of spins at each 
node into a single spin s = (d — l)/2 (d is the dimension of the corresponding space), 
as indicated in Fig. I0.4f b). That is, we apply P : Ho ® Hd — > where is a 
Hilbert space of dimension x. It can be easily shown, that this map may only reduce 
the entanglement entropy and thus we will have Ea < 2 log2(I?), an area law. We can 
do exactly the same in any spatial dimension. For instance, in a square lattice, we just 
start with four spins at each node and use a map P : H'^'^ — > H^- Thus, the states are 
generated by projecting auxiliary entangled pairs onto physical spins. The states so 
produced arc called projected entangled-pair states (PEPS) (|45|) and obviously fulfill 
the area law for fixed D. PEPS are completely characterized by the maps P we have 
to apply to each node, since the auxiliary states are in maximally entangled states 



Thus, the number of parameters required to specify a PEPS is NdD^ , where z is the 
coordination number. 

Let us consider a simple example. Let us consider D ~ 2, one spatial dimension, 
and apply the map P = |0)(0, 0| + |1)(1, 1| on each site. The resulting state is a GHZ 
state (|UT2l) of N sites. 



D 




(0.43) 



n=l 
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(a) 




(b) m IP, IP, IP, IP, IP, m IP, IP, m„ 

Fig. 0.4 (a) Dimerized state: in each node we iiave two spins wliicli are maximally entangled 
with the neighboring spins. This is represented by a line; (b) Matrix product state construc- 
tion. Starting from the dimerized state, we map the state of two spins in each site to a single 
spin. The maps P completely characterizes the state. 

If we want to approximate in terms of PEPS the ground state of Hamiltonians 
with short-range interaction, as the ones considered in the previous section, we may 
have to increase D with TV. However, since we fulfill the area law (and may only have 
logarithmic violations), this increase is expected to be polynomical with N ^ rendering 
an efficient description of the ground states In fact, Hastings has proven (|T7)) 

that in any spatial dimension and for any finite temperature, PEPS can be used to 
efficiently approximate the Gibbs state as long as the interactions are short range. 

PEPS are not restricted to spin systems or pure states. In fact, one can define 
fermionic PEPS pS)) . or mixed PEPS ([^ [50)) . and one can extend them even to 
the continuous limit ([5T|) or to the case where the dimension D ^ oo ([5^. Apart 
from that, different algorithms have been proposed and implemented in order to solve 
particular Hamiltonians using those families of states variationally [see, eg, the review 
©]• 

0.5.2 1 dimension: matrix product states 

In ID the above construction gives rise to the so-called matrix product states (MPS) 
((551 [S3]). By writing the map P„ = ^ A^^ ^ [n] \i) {a, (3\ we obtain 

d 

|*)jv= tr(^"i[l]A"n2]...A"«[7V]) |ni,n2,...,nAr). (0.44) 

ni ,. . . ,njY — 1 

Here, A"[Af] represents a D x D matrix. For translationally invariant systems we may 
choose all the matrix independent on the sites, ie 74"[Af] ~ A" [A/'] = A". 

Every MPS is invariant under the exchange of A^j — )• Xm-iA^jX^j^ , where the 
X are non-singular matrices, as it can be checked by direct inspection of (|0.44p . This 
gives us the possibility of choosing a gauge, and thus impose conditions to the matrices 
A which simplify the further calculations, or which give a physical meaning. In our 
case, we can impose 

d 

^ A"[M]U"[Af] = 1 (0.45) 

Tl = l 

as a gauge condition which implements such a choice. 

Given a MPS, one can easily determine the expectation value of product of observ- 
ables. This task reduces to calculate the trace of product of matrices, something which 
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can be implemented very easily. In fact, one can show that all connected correlation 
functions decay exponentially, and the correlation length may take arbitrary values. 
One could expect that given the construction only near neighbor correlations exists. 
However, the projections may be understood as partial tclcportations ()55|) which can 
give rise to very long range correlations. 

MPS lie at the realm of the density matrix renormalization group method (DMRG) 
Ij56p . which is widely used to determine the ground state in ID spin chains. In fact, 
this method can be understood as a variational calculation with respect to the ma- 
trices ^"[M]. The success of this method lies on the fact that MPS can efficiently 
approximate the ground states of ID chains, as explained above. 

0.6 Conclusions 

In this chapter, we have analyzed several entanglement properties of many-body quan- 
tum systems. In the first two sections we have reviewed the basic concepts regarding 
entanglement of bipartite and multipartite systems, both for pure (Section 10. 2p and 
mixed states (Section l0.3|) . and related them to more standard properties (like correla- 
tions) used in quantum many-body theories. As we have argued, whereas entanglement 
for bipartite systems is well established, for multipartite ones there exist many pos- 
sibilities of defining entanglement measures. It may turn out that one state is more 
entangled according to one measure, but not according to a different one. Apart from 
that, measures of entanglement are very hard to determine in general. Nevertheless, 
we have highlighted two measures of entanglement/correlations which are simple to 
determine and have a clear physical meaning. The first one, the entropy of entangle- 
ment, applies to bipartite pure states and is given by the von Neumann entropy of the 
reduced density operator of one of the subsystems. The second one, the quantum mu- 
tual information, measures the correlations between the two subsystems and for pure 
states it reduces to the entropy of entanglement. One can apply these two measures 
to multi spin systems in lattices by separating all the spins into two disjoint regions 
which are then considered as a bipartite system. In the particular case where we deal 
with ground or thermal states of short-range interaction Hamiltonians, the application 
of those measurements give rise to area laws, which have been the subject of Section 
10.41 These laws state that the quantum mutual information between a region A and 
its complementary one scales with the number of spins at the boundary of A (and not, 
as one would expect, with the total number of spins in A). We have also mentioned 
that possible violations to this law appear for certain critical systems. The violation is, 
however, very mild and display universal properties (independent of the lattice con- 
stant, as expected for scale-invariant states). This is very peculiar, and provides us 
with a signature of the many-body quantum states that appear in Nature in thermal 
equilibrium. In fact, the area law can guide us to find efficient descriptions of thermal 
equilibrium states, where the number of parameters does not grow exponentially with 
the volume of the system. In Section TO. 5. II we have made use of this fact and introduced 
certain tensor network states, the projected entangled-pair states, and explained some 
of its properties. In particular, in ID they reduced to the matrix product states, that 
play a very important role in a numerical methods widely used in condensed matter 
physics. 
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